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Abstract 

Using cumulant expansion for an averaged Wilson loop we derive an action of the gluodynamics 
string in the form of a series in powers of the correlation length of the vacuum. In the lowest orders 
it contains the Nambu-Goto term and the rigidity term with the coupling constants computed 
from the bilocal correlator of gluonic fields. Some higher derivative corrections are calculated. 
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1. Introduction 



One of the most exciting questions of the modern quantum field theory is a possible relationship 
between the large-distance behaviour of the confining phase of gluodynamics and string theory 
(for a review see for example^'^). The aim of this letter is to derive an effective action for the 
gluodynamics string using the averaged Wilson loop expressed in terms of field correlators^'^'^'^. 
In what follows we keep for simplicity only the lowest - bilocal correlators, which are believed to 
be dominant according to lattice data^, and briefly discuss the effect of higher correlators. The 
bilocal correlator may be parametrized in the following way^'"''^'^: 
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Here Tg is the correlation length of the vacuum, which is small in the comparison with the size r 
of a Wilson loop in the confining regime^. 

In order to derive an action for the gluodynamics string, induced by the nonperturbative back- 
ground fields, we shall consider the expression for an averaged Wilson loop using the nonabelian 
Stokes theorem^ 



< W{C) >= tr expi-^ j da^,{x) J d<Jxp{x') < F^,{x)^x, x')Fxp{x')^x' , x)>\. (2) 
^ s s ' 

At this point one should make clear of the notion of the surface S entering the integral in (2). In 
gluodynamics the Wilson loop average depends on the contour C and should not depend on the 
shape of the surface 5* bounded by this contour, and therefore one usually considers S to be the 
surface of minimal area. In what follows we generalize this definition, taking S to be any surface, 
since our aim is to derive an effective string action for the case when the surface, swept by this 
string, is arbitrary. Such an effective action contains all geometrical characteristics of the surface 
in question, and one can always specify the surface S in the final expression. 

To derive this action we shall expand the integral in (2) in powers of ^. Such an expansion 
yields in the lowest, second, order in Tg the usual Nambu-Goto term with the string tension 
proportional to the surface integral of the function D in agreement with^, while in the order 
there arises the rigidity term*'^ with the inverse bare coupling constant proportional to the first 
moment of the function 2Di — D. The sign of this coupling constant is connected with the type of 
dual superconductor, describing the nonperturbative gluodynamics vacuum^°, and hence we may 
quote the following result of the next section: if / d!^ zz"^ {2Di{z'^) — D{z'^)) < than this is a type-II 
dual superconductor (in the Abelian case the stabihty of the classical Abrikosov- Nielsen- Olesen 
strings is ensured only in the case when this type of superconductor is realized^*^). 

The main results of the letter and possible future developments are discussed in the Conclusion. 

In the Appendix, as an example of geometrical structures arising in higher orders in we 
present some of the higher derivative terms in the order Tg. The corresponding inverse bare 
coupling constants are proportional to the second moment of the function D + Di. Therefore the 
effective action of the gluodynamics string, generated by the nonperturbative confining background 
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fields, has the form of a series in powers of -f, corresponding to the series in powers of the scalar 
curvature of the manifold and more complicated geometrical structures. 

2. An action up to the order Tg 

Let us rewrite the correlator (1) in the form, which is more convenient for the future calcula- 
tions: 
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< F^^{x)^x,x')Fxp{x')^{x' ,x) >= — \ {5^xSup - 5^p5u\) 
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where D'^ denotes the derivative of the function Di by the argument. 
Our first goal is to compute the integral 



J= j d(jp,^{x) j daxpix'){5p,xS^p - 5f,pS^x) 
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Here da^,^{x) — y/g{^tfj,,^{^)d'^^, — -^s"'^{daXp){d})Xi) is the extrinsic curvature of the string 
world sheet, —2, gab — (daXp^ldbXp) is the induced metric tensor, g — det \\ gab \\,da = 



= ■^;a,b= 1,2. Expanding ^Jg{^'), txp{C), x' - x and D| 



+ Di\ 



{x-x'f 



where x' = x{^'), systematically in powers of passing from the ordinary derivatives to the 
covariant ones (without torsion) via the familiar Gauss- Weingarten formulae DaDbX^ = dadbX^— 



'^ab^cX/i 



■riinjl 



Sij, riip^daXp, = 0; i^j — 1,2, where Tab ^ Christoffel symbol, Kl^ 
is the second fundamental form of the manifold, are the unit normals to the sheet, one gets in 
the conformal gauge gab — e'^Sab 
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4Mo-T^^g^\dat^.){dbt^.) 
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where Mq = Jd^z{D{z^) + Di{z^)), Mi = J d'^zz'^{D{z'^) + D{z^)), < >= 

= tr < F^,^(0)F^^(0) > . Here we omitted the full derivative terms of the form / d^$,y/gR, where 

i? is a scalar curvature of the manifold, and used the formula {DaD'^Xp){DbD^Xp) = 

= g°'^{dat^y){dbtp_u)- The estimate for the neglected terms may be easily obtained if one assumes 

that the string world sheet is not much crumpled, so that the induced metric is a smooth function, 

which means that the typical values of ^ are of the order of r. 

Using the relations tp,xKx = 9'^{daX^){dbX,), {g^'^'g'^ + + rt'')(5aV)(^6^M)- 

■{dcXx){ddXu) = K^a^b^ — R and omitting the full derivative terms one can in analogous way 
compute the integral 

^ j da^,^{x) j daxp{x')[{x - x')^{x - x')xSyp - {x - x')^{x - x')p8yx+ 
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which occurs to be equal to 

r(i) 



+ 0{^, ), (4) 



where M^^^ =Jd^zD^{z^), Aff ^ = J d^zz^Diiz"^). 

During the derivation of the formulae (3) and (4) we exploited the fact that for any odd 

Combining together (3) and (4) we finally obtain the effective action of the gluodynamics string, 
induced by the nonperturbative background fields, in the approximation when all the correlators 
higher than bilocal are neglected: 



'^biloc. 

where 
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(5) 



a = 4T^ J d^zD{z^) (6) 
is a string tension (which agrees with^) and 

^^-T^ ( d'zz\2D,{z')-D{z^)) (7) 

is an inverse bare coupling constant of the rigidity term. 

Hence we proved the statement, announced in the Introduction, namely when 
/ dPzz^{2Di{z^) — D{z'^)) < 0, the nonperturbative Euclidean gluodynamics vacuum may be 
considered as a type-II dual superconductor, which in the Abelian Higgs Model case implies that 
the Londons' penetration depth of magnetic field is larger than the correlation radius of the Higgs 
field fluctuations, and the classical Abrikosov-Nielsen-Olesen strings are stable^''. In particular 
one concludes that when D > 2Di everywhere, the confining regime of an averaged Wilson loop 
is realized according to the dual Meissner effect mechanism^^ with the string tension given by the 
formula (6). 

In the higher orders in — more and more complicated geometrical structures, containing higher 
covariant derivatives, arise in the string action. The inverse bare coupling constants of these terms 
are linear combinations of higher moments of the functions D and Di . One can therefore establish 
some type of correspondence between the expansion of the string action in powers of -f- and a 
multipole expansion in two-dimensional gravity. A generic n-th (n > 2) term of the string action 
is proportional to some linear combination of the (n — l)-th moments of the functions D and Di, 
which are of the order of T^" < > and to the {2n — 2)-th derivative of the induced metric. 

Integrating over d^^ one obtains an estimate < > r'^(y-)^". As an example we present some 
geometrical structures, arising in the string action in the order T^, in the Appendix. 

The same effect will be due to the higher correlators. However there is one more parameter in 
the problem, which is assumed to be much less than 1 in the confining regime of the Wilson loop. 
This is the parameter of cumulant expansion^, g < >2 T^, whose n-th power estimates the 
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upper limit of the n-th order cumulant g"- < F^jjyj(a;i)$(a;i,a;2)F^2t^2(^2)--- Fnr,u„{xn)^{xn,xi) >. 
One may conclude that the higher terms of cumulant expansion will contribute to a and cuo also (as 
well as to the coupling constants of the terms arising in higher orders of — ), but their contributions 
will be suppressed by the additional powers of the parameter of cumulant expansion in comparison 
with the formulae (6) and (7) derived from the bilocal correlator. 

Therefore we see that the Wilson loop average, written through the field correlators, while 
expanded in powers of -f- gives rise to the expansion of the string effective action, which may be 
called curvature expansion. 

Note that the formulae (3)-(7) (as well as (A.1)-(A.4) presented in the Appendix) were de- 
rived in the conformal gauge, while wc dealt with the open string, sweeping the area inside the 
Wilson loop. It is known^ that in the case of a unit disc (onto which the string world sheet 
in our case may be unambiguously mapped) this gauge is accessible, but one should take into 
account diffeomorphisms reparametrizing the boundary. This reparametrization, defined modulo 
S'L(2,R) transformations, is determined by the original metric gabiC)^ ^ind thus if we quantize 
the Nambu-Goto term using the method suggested in^^, in the functional integral over all the 
metrics one should take into account not only (^-integration, but also integration over all possible 
reparametrizations. 

3. Conclusion 

In this letter we used the Wilson loop average, written via the field correlators, to derive the 
effective action for the gluodynamics string in the approximation when all the correlators higher 
than bilocal are neglected. Such an action has the form of a series in powers of the vacuum 
correlation length and in the lowest orders is given by the formula (5), containing the Nambu- 
Goto term with the string tension, defined via (6), and the rigidity term with the inverse bare 
coupling constant (7) proportional to the linear combination of the first moments of the functions 
parametrizing the bilocal correlator. It is shown that in agreement with the 't Hooft-Mandelstam 
mechanism of confinement the confining regime of the Wilson loop takes place in the case of type-II 
dual superconductor model of the gluodynamics vacuum, which in the Abelian case corresponds 
to the situation when stable Abrikosov-Niclscn-Olesen strings exist. Therefore the criterion of 
distinguishing of the confining and dcconfining regimes (or the types of superconductor in the 
Abelian case), following from the bilocal correlator, is estabhshed. In general expansions of the 
Wilson loop average, expressed in terms of field correlators, in powers of the vacuum correlation 
length and in powers of the parameter of cumulant expansion produce the curvature expansion 
of the effective action of the gluodynamics string generated by the nonperturbative background 
fields. 

However the approach suggested for derivation of this action leaves us with the conformal 
anomaly of the Nambu-Goto term if D ^ 26. One of the possible methods of its cancellation in 
D = 4 was suggested in^^ for the case of the strings in the Abelian Higgs Model. In the frame- 
work of this method passing from the field variables to the collective string ones and computing 
the Jackobian, corresponding to such transformation, one gets from it the Polchinski-Strominger 
term^^ in the action, which exactly cancels the conformal anomaly in D = 4. It should be em- 
phasized that the Polchinski-Strominger terms do not arise during our derivation of the effective 
action (5), and hence one should think about the mechanisms of cancellation of the conformal 
anomaly in D = A possibly similar to one suggested in^^. While in this letter we derived the effec- 
tive action of the string, induced only by the nonperturbative confining fields, it seems natural to 
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try to disentangle this problem in a rather elegant way taking into account perturbative gluons' 
contributions and reformulating the summation over the surfaces in the functional integraP'^^ in 
terms of the perturbative theory in the nonperturbative gluodynamics background^. Within this 
approach the perturbative gluons and ghosts, propagating inside the Wilson loop, generate the 
string world sheet excitations. The investigation of this problem will be the topic of another 
publication. 

One more set of questions is connected with the rigidity term in (5). It is known^'^, that in the 
generic case the rigid string has a crumpled world sheet, and its spectrum contains bosonic tachyon 
- the particle with imaginary mass. However as was discussed in^'^, these problems disappear if 
the /9-function has a zero at some value of the coupling constant, which may take place for example 
due to some 6'-terms in the action. This problem will be also treated elsewhere. 
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Appendix. Some examples of geometrical structures arising in the 
action in the order Tg. 

In this Appendix we present a part of the string action terms in the next order in Tg. Namely 
we demonstrate the geometrical structures arising from J in the order T^. 

One can prove that O -^—2 — in (3) equals to 



J d'Cy/^d'X Ji + J2 + J3 + J4 + P(A2)A"^'='^ 



^{dadbdcddy/g) - ]^{dat^u){dbtf,y){dcdd^/g)+ 
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■{dddedfx^) + {dbx^){dcdddedfx^)) + (dadb^/g) (^{dcddXu,){dedfX^) + '^{dcXi_,){dddedfXij,)j + 

2 f 1 1 \ 

+^{dadbdc^/g){ddX^){dedfXf,) - ^/g{datf,u){dbtf,u){^{dcddXx){dedfXx) + -{dcXx){dddedfXx)\- 



-{datnu) {dbtij.t.) {dc^) {ddXx) (dedfXx) + ^^t^^{dadbdct^^) {ddXx) {dedfXx) 



_|_ pff ^^2-^ ^abcdefij _ 



{da\^){dbX^,){dcddXi^) ( {dedfX^){didjX„) + ^{deX^){dfdidjX^)\ + {^a^b^/g){^cX^,){^d^eX^,){^fX^)■ 



■{didjX^) - ^^/g{^atf,u){^bt^,^.){^cXx){^d^eXx)i^fXp){^i^jXp) 



• {dbx,) {d^ddx,) {dex,) (dfdiX,) (djxx) (dkdixx) + I ^P'^^iX^W^^f'^^'^^^idaX,) {dbd^x,) {ddX,) ■ 
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• (dedfx^) (diXx) (djdkXx) {diXp) {dmdnXp) | . 

Here P = D + Di, \^ = ^-a^A^A*, A*i -*" = A*^..A*", and Ji, J2, J3, J4 are the following terms, which 
do not contain derivatives of P higher than of the third order and do not depend explicitly on 
\/gttiu^ \/gdattiv and derivatives of 

Ji = \^P\\^)K''''=^^{^-{dad,X^){d,ddd,dfX^^^^^ , 
J2 ^ 2^P"{X')A'^'^'^^'^^^^{dad,X,){d,ddX,) + l{daX,){dMx,)y 

■^{dedfX^){didjX^) + ^{deX^){dfdidjX^)^, 
Jz = ^V^P"(A^)A"^^^(a,x^)(a5a,.x^)(2(adaeX,)(a^a,a,x,) + {dax,){d,djd,d,x,)), 

J4 = 2v^p'"(A2)A«''^'^^^^^-^^(aax^)(a69ex^)(a,x.)(aea/x,)- 

• (^{didjxx){dkdixx) + ^{diXx){djdkdiXx)^ ■ 

Using the standard rules of computation of higher covariant derivatives and the Gauss- Weingarten 
formulae it is possible to calculate the integrals / (fXJi, / dFXJ^. The results have the form: 



{DaD^x^){2D^D^ + D'^D^ + 2DbD^D^)x^ + {D^DaX^){D^D'' + 20^0'' D'')x ^+ 



+ {DaDkD''x^)D^D^D''x^ + {DaDi,D^Xp){D''D^D^ + D^D^D^ + D^DW + D^D^D" + D^D^D^^ 

+ -{DaDbX^){D'^D''D^ + D'^D^D" + D^'D^D^ + D^D^D" + D^D^D"^ + D^D''D'^+ 
6 nc I nfe n no nc 1 n no n6 nc 1 n nft na nc 1 n na nc 7->6 1 n nft nc na\ 



+D^D''D^)x 



^D^DcDW^ + DWcD^D" + DcD^D^D" + DcD^D^D" + DcD^D^D" + DcDW''D'')x^+ 

+ -{D'^x^){DaD'^D^ + DaDbD^D^ + D^D^)x^ + -{daX^){D''D^D^ + D^D'^D'^ + D^D^D''+ 
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+D^DbD''D^ + DbD^D^'D^ + D^'DbD^D'' + DbD^'D^D^ + DbD^D^D'' + DbD^D'^D^^ 
+D''DbD^D''D''+DbD''D^D''D''+DbD^D''D^D''+DbD^D^D''D^+DbD^D^D''D''+DbD^D''D''D'')x^+ 



{D^x^)l^D^D'' + 2DbD^D' + ^D^D^ + (D^D^x^ + D^D^x^ + D.DbD'x^) 



•(3D"D^ + D^D'')x^ + {DbD^x^) i2D^D''D'' + 2D''D''D^ + -D'^D^D" + -D^D^D" + -D''D^D''+ 

\ 2 2 2 

7 \ 1 

+-D''D''D^ x^ + —{dbXf,){7D^D''D^ + ID^D^D^'D'' + ID^D^D^D" + IdD^'D^D^ + 19D''D^D^+ 
2 / 10 

+ l9D''D^D^D^+l?,D^D''D'^+l?,D^D''D^D^+l?,D^D''D^D^+D^D^D^D''+D'^D^D''+D^D'^D'')x^+ 



+{&'x^){D,D,D'^D'' + D,D^D' + D^D^)x^ - —{D,D,x^){D,D''x^){7{d^ip)d'ip + ISa^^V)- 

-^(L>„D6X^)(D"L>,x^)((aV)5V + 79^aV) - {D,D,x^){D,D'^x,) (^{d'<p)d^ip + '^d'd^cp^ - 

~{D,D,x,){D%)(^^{d^^)d'ip + 2d^d'^^ - ^{D,D,x^){D'D'^x,)(2d'^+^{d,^)'Y 

~iD,D,x,){D'^D%)(2d'^+^{d,^y^ - ^(D%)(D%)(2aV + 5(a,(^)^) + 
/or qq oq 
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(A.1) 
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J dFXJ^ = ^ 1^ {D,D,x^){D^x^){{D^x,){D''D'> + L>''D»)x, + 2(D"L>,x, + D,D''x,)(D'D' 

+D''D^)x,) + (DaDhX^) (DcDdX^) {2{D''D^x,)D'^D''x, + 2 {D''D''x,)D'^D^x, + {D'^D^x, + D^D'^x,) ■ 
.pc^d ^ D'^D'')x^ + {D''D''x^)D^D'^x^ + {D''D''x^)D'^D^x^ + (D^D''x^)D'^D''x^+ 
+(D'D''x,)D''D''x,) + {DaD,x^ + D,DaX^)(D''D,x^ + D,D%){D^x,)D''D''x,+ 

+ {DaD,x^){D,D,x, + + D'^D'x^){D'^D^ + D^D")a;,+ 

+{D'D'x, + D'^D^x,){D'^D'' + D''D-)x^) + ]^{D^x^){D^x^){DaD^x,)- 
■{D'^D' + D'D'')x, + {DaD^x^){D''D'x^){D,DdX,) (]-D''D'' + D-^DMx, + ^(D„DftX^)(D^L>%)- 



■{D,DdX,){D''D''x, 



+D'D'^x,)r-d'<p+\{d^ipf\-\^g 



^V^(D2a;^)(D%)(^2aV+ y(5c^)') - ^{DaD,x^){D- D'x^+ 



6 ric. 



+6{D,D''x^)D'D''x^ + 2(L>«DeX^)L)''D% 



6 nc„ 



1 

91 



(A.2) 



/ 



+D''D'D'')x^ + {D^DbX^, + DbD^x^) {AD'D' ^AD^D^D" + D''D')x^ + (D^x^) (4D^L>" +4Dc^"i^'+ 



+ {D^x^){D''D' + D'D'')x^] - 2(dcip)\DaD,x^ + D,DaX^)(D''D' + D'D'')x^ + 2(D^x^)D^x^]^, 



I "^'^-^^ = ^g[l{da^M^)[l3{D''D%)D% + n{D''D,x^)DW% + 7{D'^D,x^)D''D^x^ 



(A.3) 



+ 



+7{D,D^x,)D'D\] + ^{d,^f 



-{D'x^)D'x^ + 4{DaDbX^)D^D'x^ + 2{DaDbX^)D'D^^ 



X I, 



(A.4) 



where M2 = / (^2^(^2)2(^(^2) + Di{z^)), = DiD\ = did\ = did'dj&, 

{da<pr = (a„(^)(av), {dM' = a„a,9«av. 
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